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We study the distribution of equilibrium avalanches (shocks) in Ising spin glasses which occur 
at zero temperature upon small changes in the magnetic field. For the infinite-range Sherrington- 
Kirkpatrick model we present a detailed derivation of the density p{AM) of the magnetization 
jumps AM. It is obtained by introducing a multi-component generalization of the Parisi-Duplantier 
equation, which allows us to compute all cumulants of the magnetization. We find that p{AM) ~ 
AAf ~^ with an avalanche exponent r = 1 for the SK model, originating from the marginal stability 
(criticality) of the model. It holds for jumps of size 1 ^ AM < N^^'^ being provoked by changes of 
the external field by 5H = 0{N~^^'^) where A'^ is the total number of spins. Our general formula 
also suggests that the density of overlap q between initial and final state in an avalanche is p{q) ~ 
1/(1 — q). These results show interesting similarities with numerical simulations for the out-of- 
equilibrium dynamics of the SK model. For finite-range models, using droplet arguments, we obtain 
the prediction r = {d{ + 0)/dm where d{, dm and 6 are the fractal dimension, magnetization exponent 
and energy exponent of a droplet, respectively. This formula is expected to apply to other glassy 
disordered systems, such as the random-field model and pinned interfaces. We make suggestions 
for further numerical investigations, as well as experimental studies of the Barkhausen noise in spin 
glasses. 



I. INTRODUCTION 

The low-temperature response of disordered systems 
often proceeds in jumps and avalanches. I^^^^l These pro- 
cesses are beyond standard thermodynamic calcula- 
tions and are therefore relatively difficult to access and 
describe analyticalljli^Hls] jjj g, recent articl^^^, we 
succeeded in calculating the statistics of equilibrium 
avalanches (also called shocks) in a variety of disordered 
systems described by mean-field theory based on Parisi 
replica symmetry breaking. This encompasses in partic- 
ular the canonical Sherrington Kirkpatrick (SK) model 
for the Ising spin giass^^ and elastic manifolds in the 
limit of a large number of transverse dimensions. Al- 
though it has been known for a while that the equilibrium 
magnetization curve M{H) of the SK model undergoes a 
sequence of small jumps as H is increasecP^l, their statis- 
tics had not been obtained previously. The aim of the 
present article is to provide a detailed derivation of the 
distribution of avalanche sizes for the SK model. We in- 
troduce replica techniques that significantly extend the 
formalism developped in Ref . HHI to study velocity corre- 
lations in high-dimensional Burger's turbulence. It also 
generalizes previous studies of the variance of equilib- 
rium jumps to their full distribution^^BMI, expect 
this technique to be useful in several other contexts as 
well. In particular, it should be helpful to describe the 
response of complex systems to a small change of param- 
eters, a problem that arises in a variety of fields ranging 
from condensed-matter physics of complex systems, op- 
timization problems to econophysics^^HIZl 

The main result of our calculation is that the distri- 
bution of jumps takes a scale- free form, described by 
a power law of the jump size. This is intimately tied 
to the criticality of the spin-glass phase of the models 



analyzecpslj and we conjecture that such a criticality is a 
feature of a large variety of frustrated glassy systems.!^ 
The exact result obtained in the SK model finds a natu- 
ral interpretation which allows for an extension to finite 
dimensions via droplet scaling arguments. Those relate 
the equilibrium-avalanche exponent to critical properties 
of droplet excitations. 



Our results complement previous numerical simula- 
tions by Pazmandi et al^^ on out-of-equilibrium hys- 
teresis at T = in the SK model, which exhibit sur- 
prising similarities, as we will discuss. Understanding 
the relations between these results requires further nu- 
merical investigations of dynamic and static avalanches, 
both in mean-field and finite-dimensional spin glasses. 
Our results suggest to look for power-law distributed 
Barkhausen-type noise in spin and electron glass experi- 
ments, as will be discussed. 



This paper is organized as follows: In Sec. II we re- 
visit the Parisi saddle-point equations in the presence of 
a small varying external magnetic field. In Sec. Ill, we 
generalize the Paris-Duplantier equations to compute the 
moments of the magnetizations in different fields. From 
that calculation we extract the distribution of equilib- 
rium jumps in Sec. IV. In Sec. V we consider the case of 
finite-dimensional spin glasses, and using droplet argu- 
ments we obtain a power-law distribution of equilibrium 
avalanches. In Sec. VI we discuss the connection with 
previous numerical studies on spin and electron glasses, 
and propose experimental and numerical investigations. 
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II. MODEL AND METHOD 



B. Saddle-point equations 



Model and aim of the calculation 



We study the SK spin-glass model of energy 



N 

■H - - ^ Jy CrV^ 



N 



i=l 



(1) 



where the are i.i.d. centered Gaussian random vari- 
ables of variance ,P /N , that couple all N Ising spins, and 
i/cxt is the external field. 

Our aim is to follow the equilibrium state as a func- 
tion of the applied field H^^t at low temperature (3^^ — 
k^T = T <^ J. We consider small variations of the ap- 
plied field around a reference value H, i?ext = H + 

We are interested in the total magnetization in a given 
sample, 



(2) 



where F — —fcBT'liiTrexp(— /JH) is the free energy. 
Since upon variation of h of order one we expect jumps 
of the total magnetization of order \/N we define: 



1 



N 



M\ H 



h 
~N 



-dhFih) 



(3) 



where from now on we denote F(h) the free energy in 
the external field H H — 1=. Note that mi, is the sum 

VN 

of a constant part of order \/N, TFTq = ^^^^j^, plus a 
fluctuating part mh — mo of order unity. 

To characterize the statistics of these order-one jumps 
in rUh we need to compute the following cumulants in 
different physical fields hi, i — 1, . . . ,p: 

ruh^ . ..ruh^" = dh^ . . . dh^S'-^'^ {hi, h2, . . . , hp) . (4) 

It is obtained from the cumulants of the sample-to- 
sample fluctuations of the free energy, 

S^''\hi,h2, ...,hp) = {-lYF{hi)...F{hpf\ (5) 

where we denote by . . the average over disorder and 

. . . ''^ its connected average. 
These can be obtained from the generating function 
= W[h] of a — 1, . . . , n replica submitted to 
different fields /i", 



exp '■— exp 



n 
a=l 



(6) 



Note that fields ft." with replica index a are denoted with 
upper index to distinguish it from the physical field hi 
with lower index. Hence 



W[h] 



oo 



(7) 



We now derive a formula for W[h] from the saddle-point 
equations in the large- limit. 



One has: 



= ^exp /3^a^J,,a^-f/3^(i/ 



X exp 



X exp 



a ^ 



h" 

N 



Qab 



Qabcrlal 



(8) 



Note that on spins tr^, we put the replica- index a at the 
bottom, and the site index i at the top. Now we define 
the local partition sum 



(9) 



^exp 



/3''j^J2Q-b<Jaab + Pj2{H 

a^b a 



in terms of which we can write 



JV[h] 



dQab exp 

a^b 



nN 



N 
'2 



a^b 



(10) 



In the limit of iV — > oo we can perform a saddle-point 
evaluation. For /i" = this is the usual SK saddle-point 
equation in presence of a field Fl. In the low-temperature 
phase considered here, it has a set of solutions, denoted 
9a& = 97r-i(a)7r-i(b)- They are obtained from the standard 
Parisi solution qab by applying a permutation tt e S{n) 
of the indices. Each saddle point qab of the path integral 
over Qab satisfies the self-consistent equation for a ^ b: 



{^a(^b) A(q,0) - 'iab , 



(11) 



where 



refers to an average with action A from 



Eq. ([9|). Since changes in the external fields are of size 
h°'/^/N, they alter the solution of the saddle-point equa- 
tion from q = Qo to qh = qo + Hence we can 
compute the contribution to W[h] of each saddle point 
in perturbation theory. For a given saddle point, each 
contribution to e^^^^ is of the form e^t"""''!, with 



V[q,h] ■= nN 



N 
~2 



P^J^Y.'&b + NA[q,h] 

a^b 



(12) 



The saddle-point condition satisfied for any h reads 

d,^,V[qh,h]=0 . (13) 
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Using this equation, we obtain the following expansion 
in replica fields h°': 

a 

+ \Y.^^h\dU,V)[qM + 0{^) 

ab 

= V[q,0]+l3Y,h''m^ 

a 



This works, since the terms in ( f 6 1 with q < p vanish 



In the first line we used the condition (131, its total 
derivative w.r.t. and that dhqh = *^(^^) 
nate the cross term dqdhV; in the second line we used 



Eq. (11) 



The final expression for W[h] is obtained by performing 
the sum over all saddle points q^j, = 97r-i(a)7r-i(6)j 



^W[h]-W[Q]-l3J2^ h-^nm 



(14) 



after the differentiation and the ones with q > p vanish 
in the limit n — > 0, leaving the desired term q = p. 

A 




FIG. 1. The Parisi-furiction q{x), with its two plateaus for 
X < Xm and x > Xc- Note that this gives two 5-function 
contributions to the derivative of the inverse function, ^^^^ = 
x,nS{q — Qm) + (1 ~ Xc)5{q — Qc) + smooth part. 



The prime on the permutation sum indicates that the 
sum is normalized by 1 = 1. For convenience, we 
have introduced qaa to be defined later. 

Let us define the "non-trivial" part W[h] of W[h] as 

W[h] := W[h] - W[0] - /i°m^ 

a 

a 



(15) 



III. CALCULATION OF MOMENTS 



A. Generalized flow equation 



To proceed, we decouple the ha's by a Hubbard- 
Stratonovich transformation, 



JV[h] 



E>" 



(18) 



To obtain the p-th cumulant, we need to consider 
lU[{ft.°}] for p groups of ni — ain, n2 = a^n^ . . . , rip — 
apU replica with '^^^iCti — 1. Each group is subject 
to a different physical field hi, i = 1, . . . ,p. This field 
is constant within a replica group. We remind that we 
use superscript indices h°- to denote replicas, and sub- 
script indices hi to label the replica groups. The result- 
ing Wp[h] := W{hi, ...,hp):^ WUh"}] (and likewise for 
Wp[/i]) has the cumulant expansion 



^P^^^ = E ?r"' E ■ ■ • E ■ • ^^.S^'Hhn,- • ■ hO ■ 

q ^' ii=l iq = l 

(16) 

The magnetization cumulants for p > I can then be ex- 
tracted as 

ruhi ■ ■■mh/'" = dh^ . . . dh^S'^P''{hi, ...,hp) 

lim . I. dh, . . . dhWp[h] . (17) 



where /ia are Gaussian random variables with variance 
(MaA'b)^ — P'^lab, and denotes the average over them. 



The sum over permutations in ( 18 ) is equivalent to a 



normalized sum (indicated by a prime) over assignments 
{ia\ G {1, • ■ • ,p\, describing the permutation tt: 



= h. 



(19) 



nPpP nLi 



Since the permutation preserves the number of equivalent 
replica, we have the constraint ^j.ia. — naj. With this 
notation we obtain 



e^''"'l = ( E exp^/^.„MJ) . 

(20) 

As we prove in App. O this can be rewritten as 
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IZc ULi d;/» '^(ELi (^^y^) riLi ELi eMhgj'a + vi)] ' 

I 



(21) 



valid for any n < 0, and for any set of > 0, with 
Y^^=i'^i ~ 1- This identity significantly generahzes the 
formula (D6) in Ref. US 

In the case where Qab is a hierarchical ultrametric ma- 
trix of Parisi type, parameterized by the Parisi function 
q(x) with n < X < 1, the average over fia of expression 
(21 1 can be performed extending the methods of Ref. fi51 
We recall that we use everywhere ai = I and rewrite 



^wAh] _ riLi S{Y,i atVi) g{x = n, {y,}) 
iZo riLi dy* S{J2i atVi) E^ exp(yi)]" 

We have defined 

5(x;{yJ)^e-*(-^{^'» 



\a=l \i=l 



(23) 



The auxiliary fields /ii^"* have Gaussian correlations 

A'o^Vi'^ / — P'^Vlab — 9(2;)] • For convenience, we de- 

fine qaa = g(l)- 

Generalizing the method of Ref. |33]to several groups, 
we find the flow equation for the function (j){x, {yi\) de- 
fined above: 



dcj) 
dx 







2 p 



hih 



Aq{x) / d'^cf) ^^d<j> dcj) 
Ax \dyidyj dyt dy^ 



(24) 



It must be solved with the boundary condition 

0(a; = 1; {y,}) = log (j^ ^''^ ^ Hi{y,}) . (25) 

Here and below we denote y = {yi}. 

To simplify (22 1 we first evaluate the denominator. In 
the limit of n — > one can show the general formula for 
any ai with constraint a; — 1 and n < 0: 



dfy<5(^a. 

i 

J dfy,5(^a,y,)e-(-")'"-'^(^') [1 + 0(n)] 



-(-n)i-f [l + 0(n)] 



(26) 



( 22 1 to lowest non-trivial order in n, we find 



dh,...dh^Wp[h] 

_ n j &'y5{Y,,atVi)dh^_ ■ . . dhp(j){0,y) 
{-ny-P/n.c^^ 



[l + 0{n)] 



P pea 

i=i i 



X [l + 0{n)] 



(27) 



(22) Inserting this into Eq. (17) we obtain the final formula 
for the p-th cumulant of the reduced magnetization: 

= -{-TY j d^'y<5(^a,2/,)a^,...a„XO,y). (28) 



This expression is independent of the choice of a^, as it 
must be. In case the order parameter function, q{x)^ has 
a plateau for x < Xm (as happens for the SK model in a 
magnetic field iJ 7^ 0) 



{x = 0,y) = (t>{xm: V + z/3hy/qo) , 



(29) 



where • • denotes an average over z, a unit-centered 
Gaussian variable, 



see Eq. (96) in Ref. IHJ 



dz 



2tt 



(30) 



B. TBL-shock expansion 



We now solve the flow equation ( 24 ) perturbatively in 
the nonlinear term. This generates a low-temperature 
expansion which is well suited to study shocks^^. We 
write 



(l){x,y) = 4>°{x,y) +(j)^{x,y) + 
The successive terms satisfy 



dx 



/32 q(x) 92^0 



X dyidyj 



(31) 



(32) 



with initial condition (fP{x — l,y) ^ H(y). 



(3' 



dx 



hih 



dyidyj 



dyi dy-j 



Expanding also the numerator and the exponential in with initial condition (j)^{x = l,y) = 0. 



(33) 
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The leading-order equation (32) is a linear diffusion contribution 



equation, and integrated as (for x > Xm) 



4>\x,y) = H{y + zPh^q{l) - q{x)) . (34) 
Taking into account (30), we find the contribution of (ffi 



to the magnetization cumulants 



(35) 



It is shown in appendix [X] that at T = this equals 

= diif/' w = [2q{i)r' [(-i)^ + i]r(^) ^ ^gg^ 

which is a constant independent of hi. In addition g(l) — 
1 as T — 0. For p = 2 one finds at any temperature the 



-J,c,(0) 



9(1)- 



(37) 



Even though at T = this is the full result for the sample- 
to-sample fluctuations of the magnetization, at finite T 
there will be an additional piece from (fx^ obtained be- 
low. Similarly, to obtain the full finite-T expression of 
higher-order moments of m^^ , contributions from ^p^*^ 
are needed. However, here we focus on T = 0. 

We now turn to the calculation of the contributions 
which capture the information about jumps, which are 
of order 0{\hi — hj\) in the limit T — 0. It is contained 
in the contribution of </>^ and only in that contribution, 
as was discussed in Ref. Ell Higher-order functions 4^ 
contain contributions of order 0{\hi — hj \p) at T = 0, en- 
coding information of multi-shock correlations. To first 
order in the non-linear term we find, extending the cal- 
culation in Ref. [31] 



2 ^ x' 



f /5 q(x ) DJJ / -* \ QH / ^ \ 



(38) 



As in Eq. (30), ry, zi and Z2 are independent unit-centered 



Gaussian random variables, and D^'x ■= \/ <l{x') — q{x). 

We now change integration variables from a; — > g and 
define x{q) := x{q)/T and h := h/T, the "thermal 



boundary layer variable"'2il for the external field. Using 
Eq. ( 28 1 the contribution of the first-order term to the 



magnetization cumulant becomes, denoting qc := q{xc), 
and qm-^q{0), see Fig. [l] 



mfi, ■ ■ ■ rUh^ 



T f°° ^ / \ 7 ~ \ — 7 ~ — 



■A+,A- 



• i=l 



H(y + hA+] - H(y + hA 



(39) 



A± are centered Gaussian random variables with cor- It is convenient to introduce 
relations defined from the above independent Gaussian 
variables as 



A+ ^ 7]y/q^^ + z:^^/q^^ (40) F := A+ + , G := A+ - A^ (42) 

A- = vV 1 ^ Qm + Z2\/qc — q ■ (41) in terms of which one can integrate by part 
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qc -oo -oo / pi q2 

''.-.(i) = (_l)P9-^...a-^|; f Aqx{q) f dF f dG [dl - dl f""^^ ~ 4[g.+,-2g„l - 4[,.-,l 



— oo —OO 



P poo 



n ^yAT. 



2ix^2{qc + q- 2q,n]2{qc - q) 

2 



H ( V+ + G)]-H (y+ ^fi{F - G) 



(43) 



The differential operator 9a+<9a_ — dp — Oq acts only on the Gaussian measure. Note that its action is equivalent to 
dp — dQ = d/dq. One can thus integrate by part over q to get 



rrihi . . . rrihj, 



dq 



n/ d..<E 



H ( ij+ \kF + G)\-H (y+ ^hiF - G) 



The measure over F and G is defined by 



— oo — oo 



f{F,G) 



F.G 



dF / dG 



exp 



F G 

[q,+q-2qr^] ^ 4[ge-9] 



-OO — OO 



2Ti^j2[qc + q- 2q.m]2{qc - q) 



f{F,G). 



Equivalently, one can write in terms of and A. 



with measure 



3m * 



H (y + hA+] - H [y + hA 



fiA+,A^) 



= / d^ , 



dA^ 



exp 



+q-2q^] ilqa-q] 



'\j2{qc + q- 2q,n]2{qc - q) 
I 



f{A+,A^). 



(44) 



(45) 



(46) 
(47) 



The boundary terms in the integration by part vanish, recovering the form obtained in RePi'. For T > and 
provided that whenever q{x) exhibits a plateau for < h2 hi one finds 
X < Xrn it is included as a <5 function. 



Using the expression (25) for H{y), the formula (44) 
allows us to compute the thermal boundary-layer form of 
the p-th cumulant. We give here the result for p — 2: 



-,7,c,(l) _ 



q{x) dx - q{l). 



(49) 



Added to Eq. (37), this gives the correct total sample-to- 



■J,c,(l) 



1 dx{q) 



dG 



exp 



xG^(/ii - /i2)coth 



y/ATr[qc - q] 

{hi - h2)G 
2T 



sample fluctuations of the magnetization. The fact that 
higher terms (j}^ do not contribute to this variance can be 
verified by a direct expansion of Eq. (24) in q{x). 

For general p we only study the limit T 0. 
For convenience we introduce the notation Am = 
max(A+,A_) - (^^+|G|)/2, A^ = min(A+,A_) = 
{F — |G|)/2. The calculation is performed in Appendix 
(48) [Bland we obtain for p > 2: 



rrihi . . . ruh^ 



^,=.(1) = t I ^(Am - A„) -hpAf,, + J2 {K-m+i ~ /ip-™)A?'-"ASJ + h,Al, 



-A+,A- 



1 . dxiM) , 

2./„ dq 



(50) 



Note that we have put back the physical field h = Th, making evident the result in the limit of T — > 0. 

I 



As an example, for p = 2 we obtain 
1 ri 

4 



which is the T = limit of ( 48 ) . Note that this describes 



mh^ruh, 



dq\h,-h2\'^\Gf m 

dq 
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the correction to order |/i2 — /ii | to the 2-point function of 
the magnetization ( 37 1 . This encodes the second moment 



of the jump-size distribution, as was discussed in Ref. HH 
for the random-manifold problem. We now turn to the 
determination of the full distribution from the above cu- 
mulants. 



we have to show that 



(M(/ii)---M^)-MO)---A^(o)) 



(54) 



C. Distribution of jumps 

We now derive the distribution of jumps by showing 
that the above result is identical to a p-point correlator 
of the magnetization of a two-level system, whose char- 
acteristics (jump size and jump location) are distributed 
in a simple manner. We notice that the above expression 



( 50 1 for the cumulants can be reexpressed as 



_ 1 
^ 2 



Aq 



d^(g) 



(52) 



x|G| {^i{h,)■■■^i{hp)~^i{o)■■■^lm^^ 

in terms of the "random magnetization" variable 
^x{h) = 9{h - he) Am + OiK - h)Am 
h. 



F,G 



F 

- + sign(/i 



2 



(53) 



It exhibits a jump of size \G\ at location h — uni- 
formly distributed on the real axis with unit density. F/2 
is interpreted as the mean magnetization. To prove this, 



where the overbar denotes averaging with respect to he- 
To see that, consider an interval [/iq, /il] containing and 
/iQ < < /ii < ... < /ip < /iL, in which he is uniformly 
distributed with density po = 1- 

The probability of h^ to be in the interval 



is po(^p-m+i - hp_rn)- Then 

gives the corresponding term in 
the sum ( [54[ ) (multiplied by po)- At the edge, 
when he < /ii, it gives po{hi — hQ)A^ 



p(/ii) • ■■ ^i{h, 



> hp it gives poihL 



When 



Subtracting 



M- 



p(0)---/i(0) = PohLAf^ - PohnAlj yields (j^, (mul- 
tiplied by po: which is set to unity). Both mean 
magnetization F/2 and jump-size \G\ are obtained from 
Gaussian variables with a q-dependent variance. The 
full result in (52) is obtained by weighting with the 



probability distribution of the various values of q. 

A given shock at field h = hg \s characterized by 
its magnetization jump of size Am ~ — Ti/j- 

(always > 0), and its mean magnetization rris — 



at the shock. From the above 



M{H+^) + M{H)\ 

we can extract the joint density (per unit interval of 
/i), p(Am,(5m), of shocks of size Am, and shift 5m — 
nis — rriQ. It is defined as: 



p(Am, Sm) = lim — 6\ Am — 
hio h \ 



M{H 



h N 



M{H) \ 



N 



6 Sm — 



/N' 



M{H) - 2M(iJ) j 



2\/N 



(55) 



and can be extracted from Eqs. (52 1 and (451, identifying 5m with F/2, and Am with \G\ as discussed above. This 
leads to 



p(Am, (5m) = 0(Am)Am / Aq ^^' ^ 



+q--2q„ 



j)exp( 



{Amy 

4 [go -9] 



^/TT[qe + q- 2q„,] ^47r[gc - q] 



(56) 



Note that after integration over q, the jump size and the 
magnetization shift become correlated. Integrating out 
the magnetization shift, we obtain our main result, the 
density of shock sizes per unit interval of h (at T = 0) : 



p(Am) 



e(Am)Am f'^dq^^—, 



/_ {Araf \ 



q] 



(57) 

This formula is valid for a large class of models described 
by replica symmetry breaking saddle points, as empha- 
sized irfi^. Here, we will focus on its application to the 
SK model. Apart from the prefactor Am, the above for- 



mula is essentially a superposition of Gaussians (at fixed 
overlap distance q). Their contribution is weighted by 
the density 



/ ^ dx((7) 1 



(58) 



where P{q) is the sample averaged probability distribu- 
tion of overlaps between metastable states sampled from 
the Gibbs distribution.'^ The weight v{q) can be inter- 
preted as the probability density of finding a metastable 
state at overlap within [g, q -\- dg] and energy within 
[E,E + dE], with E close to the ground state. We will 
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come back to this interpretation below. 



magnetization jump. It is 
p(Am)AmdAm — 



A useful check of Eq. (57) is provided by the average 

dx{q) 



q{x 



dq- 



dq 



-[qc - qix)] 



lim — 

T^O T 



dx\ 



(59) 



where we remind that our definition of x{q) contains a 
(5-function contribution at each plateau, so that the final 
integral over x runs again from to 1. This formula is 
generally valid^S. For the SK model it can be rewritten in 
terms of the thermodynamic (field cooled) susceptibility, 



p(Am)Ar7idATO — lim fypn 



1 {T=0) 
XZFCj — XfC ' 



since in the SK model the intra-state (zero-field cooled) 
susceptibility, xzfc = [1 ^ Qc]/T, vanishes linearly as 
T — >■ 0. Thus, the response is entirely due to inter-state 
transitions in the form of avalanches (shocks). This is 
in contrast to other mean-field models, where even at 
T = part of the response is due to smooth intra-state 
polarizabilitjff^. 




0.01 0.05 0.10 0.50 1.00 

FIG. 2. Power-law density of jump sizes for the SK model. 
The power law receives contributions from all overlaps 1 — q. 
The curves in the lower part show the contributions from 
[1 — q) = 2~'' , k = 1, 12, each of which takes the form of 
the jump density in mean field glasses with 1-step replica 
symmetry breaking. The three nearly coinciding lines on 
the top show p(Am) evaluated from Eq. (57 1, for external 
fields H = 0,0.25 and 0.5, respectively, using approximations 
for q{x) described in the text. The increase of H decreases 
the cutoff at large Am, while the avalanche distribution for 
Am ^ 1 is a universal power law, not affected by H. 



IV. APPLICATION TO THE SK MODEL 

A. Study of the distribution of jumps: H = 

In order to evaluate the distribution of jumps, we need 
the full replica-symmetry breaking solution of the SK 
model in the limit of T — >■ 0. The increasing function 
q{x) is well characterizecP^lt^H^ even though no closed an- 
alytical formula is known. q(x) has a continuous part 
up to the "break point" Xc ~ 0.55, and is constant 
for Xc < a; < 1. In the limit of T this con- 
stant qc behaves as 1 — 1.592r^, and q{x) becomes es- 
sentially a function of a; = x/T that we call q{x). In 
the absence of a magnetic field H — 0, q„i — and 
q{x) w ^ with 1^(0) = 1.34523 at small x. ^8 At 
large x the function crosses over to the asymptotic be- 
havior 1 - q{x) « 4CVx2 + BT^ with C = 0.32047 and 
B = 0{1). This leads to a power-law tail for the weight 
of large overlaps q — )■ fl^, 

-3/2 



zy(g|l> l-g>r^) =.C(l-g)^ 



(60) 



We can now analyze the jump-size density using formula 
(57). We obtain analytical expressions in the limits of 



small and large Am. Numerical calculations describing 
the full range are shown in Fig. [2j For small Am the 
integral over q is controlled by 1 — g ^ 1, and we can 
approximate 

Cdq —--^ 



p{Am) 



-co (1 

2C 1 
0r (Am)^ 



q)3/: 



-Am- 



V 4(l-g); 

^47^(1 - q) 



exp 



Am < 1 



(61) 



with r = 1. The universal exponent t — 1 for jump sizes 
7V^i/2 <^ Am <C 1 results from superposed contributions 
from many overlaps, i.e. all scales, as illustrated in Fig. [2] 
The asymptotics for large Am is controlled by small 
q <^ 1, i.e., by transitions between very distant states. 
Injecting the density of states i/(0) near q = and Taylor 
expanding in q inside the exponential yields the estimate 



p(Am) w i^(0)Am / dq 
Jq 

~ 0f (Am)^' 



exp 



(Am)^(l+g) 
4 



V47r 

Am > 1, (62) 



with r' = 1. We see that avalanches with Am ^ 1 
(Ail/ ^ V N) are exponentially suppressed. 

Plots at intermediate Am = 0(1) are shown in Fig. [2] 
for three different values of the external field. As no 
analytical closed form for q{x) is available, we have used 
approximations of the type i(g) = {aq+bq^)/y^l — q with 
a = 1.28 and b = -0.64, proposed in the litcraturePESl, 
and a sharp lower cutoff q^in{H) = l.Oi?^/^. 



B. Distribution of jumps: H ^ Q 

In the presence of a finite field iJ, Parisi's solution 
develops a plateau at low x: 



q{x < Xm) = qm{H), 



(63) 



where q„i[H) « l.QH^/^ and Xm ~ ^{^)Qm{H) for small 
H, while q{x) is nearly unchanged for x > Xm^^^ '^^ It 
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is convenient to rewrite formula (57) as 



/5(Am) = 0(Am)Am 



exp 



\ 4.{ga-q)J 



dq^q) I—. r 



(64) 

where the density of states v(ci) contains a piece 6{q — 
qm)xm/T when q{x) exhibits a plateau at a: < Xm, hence 
the notation g~ in the integral. Thus, the effect of a 
magnetic field is to change the behavior of the jump dis- 
tribution at large Am, where it is now dominated by the 
plateau: 



exp I 



{Amy 



AH)] 



(65) 



Comparing with Eq. ( 62 ) we find an effective exponent 



t' = — 1 (instead of 1) in the tail of the distribution. The 
formula ( 65 ) holds only if we can neglect the contribution 



of the continuous part of q{x). A simple comparison with 
the previous section shows that this holds when Am 3> 
H~^/^. For 1 < Am < Am// the 



1/2 



behavior crosses over to a formula similar to ( 62 ) with 
t' = 1. 

Note that a small random field also produces a plateau 
in q{x), and hence, we expect its effect on p(Am) to be 
rather similar to that of a uniform field. 



C. Interpretation for the SK model 



To find a natural interpretation of formula ( 57 1 we con- 
sider what happens upon increasing h from hi to /12. If 
we take ~ h2 — hi <^ 1 we only need to consider 
the possibility that the ground state and the lowest-lying 
metastable state cross as we tune h, corrections due to 
higher excited states being of order 0{h^2)- 

We now argue that the disorder-averaged density of 
states of this two-level system is given by i'{q)dqdE, 
where i>{q) was defined in Eq. (58 1. Indeed, the defi- 



nition of the overlap distribution ^((7) is 
P{q) X! WaW~f5{q - qa^ 



(66) 



where Wa = exp(— /3Fq)/^^ exp(— /3F^) is the Gibbs 
weight of the metastable state a. At low T we can restrict 
to the two lowest states, which yield the leading-order 
term in 



Here v{q, E) is the joint probability density of overlap q 
and free-energy difference E between the ground and first 
excited state. Hence Eq. (58) holds with v{q) = j/(g,0). 

The two states differ in — N{l — q)/2 fiipped spins. 
In the SK model the magnetization is uncorrelated with 
the energy, and one thus expects the magnetization dif- 
ference between the states to be a Gaussian variable of 
zero mean and variance (at fixed overlap) 



{Am^)g = ANfi/N ^2{l-q). 



(69) 



When h increases the energy difference between the first 
excited state and the ground state changes from E (for 
h = hi) to E — ft,2iAm, where /121 := h2 — hi > 0. 
Thus, if Am > 0, a jump at equilibrium occurs when 
^21 — E/Am. For the shock probability per unit h one 
thus expects 

rqc poo 

p(Am) = lim / dq dEv{q,E) 



exp 



(Am) 
2{ArrP 



V27r(Am2), 



5\ h2i 



E 

Am 



(70) 



reproducing Eq. (57) upon integration over E. 



This argument strongly suggests that the joint density 
(per unit of h) of jumps with characteristics q and Am, 
is given by 



p(Am,q) = 6{Am)Amh'{q)- 



exp 



V 2(Am^)J 

^2^(Am2), 



(71) 



Integrating over Am we find the density of jumps with 
overlap q as 



or for the density of flipped spins, A^h = '•^ f'^ 

2Nn 



(72) 



N 



(73) 

Let us now consider avalanches with Nfi ^ N. Using 
Eq. (leOl), we find 



Piq) 

and the power law density 



C 1 

y/nl-q' 



(74) 



(75) 



P{q) = {l~T)6{q-qc)+Tpi{q) + 0{T^) (67) with p = 1. 



as 



= Tj^{q,0) + O{T'^). 



(68) 



D. Comparison with numerical work 

For the SK model, there is no numerical study of equi- 
librium avalanches to date. However, in a pioneering 
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work, out-of-equilibrium avalanches at T = were stud- 
ied numerically along the hysteresis loojP^', and found to 
exhibit criticality, i.e. a power-law distribution of mag- 
netization jumps. The external field H is increased adi- 
abatically slowly until a single spin becomes unstable. 
The latter is flipped and triggers with finite probability 
an avalanche of further spin flips, during which H is kept 
fixed. The typical difference in applied magnetic field be- 
tween adjacent jumps scales as N~^/^, which is the same 
scaling as in our calculation. During the avalanche a se- 
quential single-spin-flip update was used to ensure the 
decrease of the total energy. Interestingly they observe 
the same scaling of the jumps of total magnetization, 
AM ~ 7V^/^, and the number of spin flips (which we as- 
sume to be of the same order as the number of spins that 
have flipped an odd number of times), iVfl N, as in 
our present calculation for equilibrium. It is interesting 
to note that this implies that a typical spin flips on the 
order of N^^^ times along one branch of the hysteresis 
loop. A very similar density of avalanches with the same 
exponents t = p — I and a crossover at Am ~ 1, as 
analytically obtained for the statics here, was observed 
in the numerics. This similarity is surprising since the 
states reached along the hysteresis curve are quite far 
from the ground state, as evidenced by the width of the 
hysteresis loop. Nevertheless, the visited states share an 
important feature with the ground state: self-organized 
criticality. Indeed, the distribution of the local fields 
hi = Jij'^j + the energy cost to flip spin 

i only, is observed to display a linear pseudogap"^^ as in 
the equilibrium"^^, marginally satisfying the minimal re- 
quirement for metastability. 

To understand better the relation between static and 
dynamic avalanches in the SK model, it would be useful 
to perform both equilibrium and dynamic simulations. 
In particular, it would be interesting to determine the 
prefactor of the power-law for the density of jumps, which 
we have computed here for equilibrium, but which has 
not been determined in Ref. [511 because they normalized 
the jump density. It would also be interesting to compute 
the probability density of overlaps between states before 
and after an avalanche, and compare with the expression 



( 71 1 derived in equilibrium. 



One could measure the joint density of overlaps and 
avalanche-sizes, 



PH{Am,q) := {S[q-1 



N 



S At 



AM\ 

7W) 



H 

(76) 

where the average is taken for fixed external magnetic 
field (i.e., in practice for H E [H — 6H, H + SH], with 
SH small). It would be interesting to check whether this 
joint density takes a form as in Eq. ( [tT] ) with {Arn?)q — 
2{1 — q). In this case, this might allow to define a dy- 
namical overlap-distribution ^{q) (to be interpreted as 
the T = limit of Pdyniq)/T). 



V. DROPLET ARGUMENT IN ANY d 

Let us now discuss the Edwards- Anderson model in di- 
mension d. We first give a scaling argument to predict 
the avalanche exponent based on a droplet picture. Sub- 
sequently we will show how the previous result for the 
SK model can be recovered and interpreted in the same 
spirit. 

To determine the first avalanche as the field is in- 
creased, we need information about the lowest-energy 
excitations of a given magnetization, which will scale in- 
versely with the volume. More precisely, we expect the 
lowest excitation energy for a droplet-like excitation of 
linear size L to scale as 



1 



(77) 



This is argued as follows: Standard droplet argument^^H 
stipulate that the lowest-energy excitation of linear size 
L, including a given spin, grows typically as . These 
droplets are in general objects of fractal dimension df < 
d. We thus assume that one can cover the system of 
volume V by V/L'^' droplets, and that they are uncor- 
related. This implies the scaling (77) for the droplet of 



minimal energy. The density po of such single-droplet ex- 
citations near the ground state thus behaves as podE — 
dE/E^iniL), or Po = uo/L^ x l^/L* 

The magnetization jump associated with the overturn 
of a droplet of size L is assumed to scale as L'^"'. Of 
course, dm < d[. The numerical stud suggests that dm 
is rather close to df. We assume the total magnetization 
of droplets of size L to be uncorrelated with the energy, 
and distributed as Pl{AM) = L-'^^^m{AM / L'^^). In a 
vanishing field, low-energy droplets are believed to exist 
at all length scales. 

We make the standard assumption that droplets at 
scale L are uncorrelated from droplets at scales > 2L. 
By analogy with the reasoning given for the SK model, 
one argues that the density of avalanches per volume, per 
unit field H , and per unit magnetization change AM is 
given by 



p(AM) « lim -!- 
SHia V 



dL 



dE 



-S'min 



xSlSH - 



E 



AA/) 



Pl{AM). 



(78) 



Using the above expressions one finds 



p{AM) 



1 



dz ^Ijm{z)z'' , 



{AMY dm JO 
valid for AM ^ 1, with the avalanche exponent 

df + 



(79) 



(80) 



This prediction is very general. As discussed in Ref. [161 
it also gives reasonable predictions for elastic interfaces 
in random media. The formula was recently rediscovered 
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in the context of the ferromagnetic phase of the random- 
field Ising modeP^, in which case dm = df, and thus 
r = 1 + 0/cifP 

It is interesting to point out the close analogy be- 
tween the exact expression for the SK model ( 57 ) and the 
heuristic droplet argument ([78]). In the SK model the role 
of spatial scale is played by the overlap distance 1 — q, and 
the logarithmic sum over scales J dL/L goes over into an 
integral dq/{l — q). The equivalent of £^min(i) is given 
by the typical gap at distance 1 — q, which is known to 
b^^ Aq = (1 — g)^/^. Finally, the distribution of mag- 
netizations at fixed droplet scale, Pl(AM) is given by 

^ (An 

exp ( — ^ — 

Pl(AM) 



4[l-g]- 



Putting these elements to- 



gether and substituting them into Eq. (78) without the 
volume normalization factor, one recovers expression (|70p 



with vlq) given in (60). Note that changing variables 
from H to h = N^l^ and AM to Am = N-^/'^AM 
does not change the density of avalanche sizes per unit 
field and unit jump size. 

In the presence of a finite field iJ, droplets are believed 
to be suppressed above a scale Lh ~ l/H'' (with 7 > 
0). This implies that integration over droplet scales in 
Eq. ([78]) is cut off at Lh leading to 



(81) 



1 y roci 

p(AM) = — — dzi;M{z)z\ 
{AMY dm Jaj\//l^" 



which cuts off the power-law decay of the avalanche-size 
distribution at AM ~ L 



H 



At small but non-zero temperature we expect several 
effects. First, there is a thermal rounding of all the mag- 
netization jumps, which is apparent in Eq. (48) and was 



discussed there. The equilibrium jumps are smeared out 
over an interval Ah ^ T/y/TxFC- In order to be distin- 
guishable from the sample-averaged increase of magneti- 
zation, the avalanches should be bigger than the latter 
Am AhxFC ^ V^Xfc ^ T. Above this scale, the 
avalanche distribution is unchanged for T <ti Tc- 



VI. CONCLUSION 

We have introduced a method based on replica tech- 
niques to compute the cumulants of the equilibrium mag- 
netization in the SK model at different fields. From their 
non-analytic part we have extracted the distribution of 
magnetization jumps at T = 0. It exhibits an interest- 
ing power-law behavior, characteristic of the criticality 
of the spin-glass phase. We have also obtained a predic- 
tion of the avalanche-size exponent for spin glasses in any 
dimension using droplet arguments. We have compared 
with numerical simulations of the out-of-equilibrium dy- 



namics of the SK model and found striking similarities 
with the static calculations presented here. 

It would be very interesting to investigate avalanches 
in small fields in realistic models, as the finite-range 
Ed wards- Anderson model in 2 and 3 dimensions, to test 
some of the predictions that we obtained using droplet 
arguments. Furthermore, experimental measurements of 
power-law Barkhausen noise in sp in glasses (e.g., by mon- 
itoring magnetization bursta^ISl) could provide comple- 
mentary insight to earlier investigations of equilibrium 
noise^. 

We expect similar critical response upon slow changes 
of system parameters in many other systems described 
by continuous replica-symmetry breaking, as, e.g., in var- 
ious optimization problems (minimal vertex cove^^, col- 



oring^, and fc-satisfiability ° close to the satisfiability 
threshold, and in the UNSAT region at large k. Like- 
wise, in models of complex economic systems, one ex- 
pects a power-law distributed market response to changes 
in prices and stockJ^H. Avalanches have also been pre- 
dicted to occur in electron glasses with unscreened l/r 
interactions, and have been studied numerically in detail 
in Ref. [25] They find an avalanche exponent r = 3/2, 
which is reminiscent of the value found for disordered in- 
terfaces and random-field systems at the upper critical 
dimension. 

Finally, we comment on possible future avenues to ex- 
plore. It would be interesting to study analytically the 
dynamics of avalanches in the SK model. In principle 
one could use methods developped for the aging dynam- 
ic^^. In the simplest framework, one studies relaxation 
from a random initial state, in which case the overlap be- 
tween initial and final state vanishes at large time. Pre- 
sumably the hysteresis cycle selects a sequence of states 
which have non-trivial subsequent overlaps. This remains 
a challenge to describe analytically. A more modest, but 
still non-trivial goal consists in describing the dynamics 
starting from an equilibrium state upon an increase of 
magnetic field by a small amount ^ N^^^"^ . 

It would be interesting to study whether the states 
visited dynamically along the hysteresis curve, and the 
avalanches triggered, have a relation with the marginal 
TAP states at high energies and their distinct soft 
mode^^. It would also be interesting to analyze the 
multi-shock terms 0(|/i|'^'^^) in the magnetization cumu- 
lants, allowing to determine whether there are correla- 
tions between successive jumps. 

We thank L. Cughandolo, S. Franz, M. Goethe, 
M. Palassini, G. Zarand, and G. Zimanyi for interesting 
discussions. Wc thank KITP Santa Barbara for hospital- 
ity, while various parts of this work were accomplished. 
This research was supported by ANR grant 09-BLAN- 
0097-01/2 and in part by the National Science Founda- 
tion under Grant No. NSF PHY05-51164. 
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Appendix A: Zero'th order cumulant for the magnetization 



Here we evaluate the contribution of (fp , Eq. (35). At T = 0, one can set H[y\ — >■ maxi{2/i}, to simplify to 

-J,c,(0) _ 



"i/ii . . . 



i 

/P 
i 4=2 

z 

/P 
dPy5(^J2a,y,^ zP Y[6(yi ~ y, + Pz^q{l) - q-^[hi - h,]^ 
i i=2 
z 

/P 
dyi S(^yi + '^ ai(3zy^q{l) - q„i[hi - hi]j zP 

i=2 

= m - q„r' ^ = [mi) - ^.n)]'" '^''^'t'L'^^'^^ , (Ai) 

which is the result given in the text. 

Appendix B: Magnetization cumulants to first order in the shock expansion 



Consider formula (46). In the limit of T — > 0, h = h/T becomes very large, and we can approximate H{y) = 
maxi(yj). The idea of the following calculation is that taking a field derivative yields a derivative of H[y), which is a 
i5-function, eli mina ting one integration. 



To evaluate (46), we start with the cross-term, and choose without loss of generality hi < h2 < ■ ■ ■ < hp-. 

A+ ,A- 



1=1 

where we have denoted Am ■= max(yl+,y4_) and Am ■= min(yl+, The derivatives can be written as 



{-l)Pdi^^...df^^'^ 1'"' ^'^^n dy^s[j2^^y^)H{y + hAM)H(y + hAm) , (Bl) 



'h [y + hAM) H [y + hAm)] ^ J2 E % ■ ■ ■ + f^"^^^) ^h,, ■ ' ' ..^(^^+^^" 

'"=0 0.}, {/c,} 

(B2) 

where the sum is over partitions of the p fields hi into two groups of m and p ^ m fields with ji < ... < j,n and 
ki < . . . < kp-„i- The multiple derivative (with at least one derivative) of the first factor of H can be written as 

df^^^ . . . d-^^^ H(y + IAm) = i-ir^^AZ n ^iVji + ^ji^M - Vje - hj,AM) \{ QiVj, + Hj^Am - Vk. - KAm), 

1=1 i=l 

(B3) 

This equation is proven by noting that 

p 

1. max(yi, . . . ,yp) = ^ J| 6'(?/,; - yi). 

i=l l^i 

2. dy^ max(yi, . . . , y^) = JJ^ 9{yi ~ yi), since derivatives of the 6'-functions cancel in pairs. 

3. a further derivative of 6{yi — yi) w.r.t. yi gives —5{yi — yi). 

This result is a consequence of the fact that the maximum of m variables depends on p < m variables if and only 
if these are mutually equal. We note that this expression is symmetric in the {/ij^, . . . , /ij^}, and that a similar 
expression holds for the second factor. 
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The terms m = and m = p have to be considered separately, which we do now, starting with m = p: Using ( B3 ) 
and ehminating all the (5-functions from the derivatives of H yields 

n / ^^y^^ {j2 '^'y^ ^h,--- \H(y + KAm) H(y + 

= (-1)^-^*;,^ j dj/i 5 ai[yi + KiAm - KAm]^ max, + KiAm - kiAu - ^m)} 

where to get to the last line we have used ai — 1 and imni{hi} — hi. 
Likewise the term to = gives 



(B4) 



n S ^^y^ H{y + IAm) df^^... df^H^ 



y + hA„ 



(B5) 



Let us now discuss the terms m = 1, . . . ,p — 1. Consider 
n dy, <5 ( ^ 9^^^ . . . df^^^ H (y + IAm) df^^^ . . . df^^^_^^ H {y + hA,. 



/oo poo pin 

dj/ii / dyk, n ® (%i + ki^M ~ Vk, - ^fei^m - {Am - A„,)hk, 
-OO J — OO i—\ 

m 

^ n ® (~ + hj^AM ~ 2/fei - hk^Ayn + [Am - An)hj, 
1=1 

X 6 ^^ai?(2/ji +AM{hji - /ijj + ^ a,(yfei + A.m{hki ~ hk,)^ 

/OO pOO 
dyj, / dyk.e (y^, + hj^Au - yk, - hA-n - {Am - A„Omax,=i,...,p_„/ife^ 



X (- \yn + hj^M - Vki - /ifciA,„j + {Am ~ A,n)wxii^x„„^rahj}^ 
X S (^^ai{yj^ +AM{hj^ ~ hj^) + ^a,{yki + Am{hki -^-fc.)^ 



(B6) 



Note that by going from the first to the second line, we have used the (5-functions to fix j/j, = yj-^ + (/ij^ — hj^)AM, 
and yki — yki + {hki — ^fei)^m- From the second to the third line we have used that Am — Am > to simplify the 
products of 0-functions. 

The product of the two O functions implies that the contribution is non-zero only if the partitions satisfy hji. > hk^ 
for all z, £. Since we ordered hi < . . . < hp, this identifies the set of hk^ to be {/ii, . . . , hp-m}, and the set of hj^ to be 

{^p— m+l ; • ■ ■ 1 ^jj} - 



Making in (B6| the shift of variables yj-^ — > yj-^^ + yk^ eliminates yk^ from the Q functions, and allows to do the 
integral over the latter, resulting into 

([B6f = (-1)p-2a™ ^^7™ f dyj.Q + h,,AM - hk^A^ - {Am - A^) max,^i,...,p_™{/ifcj) 

J — OO 



xe 



Vji + hjiAn - hkiA„ 



+ {Am - Am) mmt=i^...,m{hji} 



{-ly-'^AZA^-^ [ Ay, My,, - {Am - Am)hp-^-i) + {Am - Am)hp-rn) 

J — OO 

{^l)PA'j^jAPm'''{AM - Am) (hp-m+l " hp^m) 



(B7) 
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Putting all terms together, (Bl) becomes 

P roo 



i-^r n dy. '5( ^ aW.)dk^ ■ ■ ■ df^^ [h {y + Um) H {y + Ia. 



p-i 



= -{Alt' + ^'rrt')h + (Am - A™) -hpAP^ + J2 Chp-m+1 - h,_^)AP^"^AZ + h.A^M 



(B8) 



where h := '^^aihi. The first term ~ h disappears once we subtract the contributions from the non-crossed terms 
l[H{y + hAM)H{y + hAM)] and l[H{y + hAm)H{y + hAm)]- This leads to the formula ( 50 1 given in the text. 



Appendix C: Proof of Eq. (21) 



Here we prove that for all sets of Ha with replica indices a = 1, n the identity 

IZo nr=i dy» (^CELi riLi ELi eMh^^J■a + vt 



X! exp [ ^hi^fia ] = 



(CI) 



holds. By definition of the primed sum, the left hand side reduces to 1 for fia = 0, in which case the idenitity is trivial. 
We now prove the identity by series expansion in Ha ^ 0. 



We define 



K^ifia) := -J 



exp (hifia) 



^ Ej=i exp (hjfia) 



1, 



(C2) 



which has the property that Ki{fia = 0) = 0, as well as 
J2^=i Ki{iJLa) = 0. We can then write 



exp{hiHa) = [1 + Ki{fia)] - ^ exp (hjHa) 



(C3) 



Analogously we define 



1 



■^{y} — -y'exp(2/,), 

rj ^ — ^ 



i=l 

exp(y») 

my) 



(C4) 
(C5) 



so that = 0, and exp(yO = Af{y) [1 + A,{y)] . 

With this one finds 



(C6) 



j^^y) [1 + My)] [1 + KM] - E {h.fia) 



P 



1 ^ 

l + -^A,(y)if,(/ia) 



^ i=l 

p-1 



1 + E 



A,(y) - Ap(y) 



i=l 



r 



With this notation the identity (CI ) to be proven can be 
restated as 



E 

.e{i,...,p}|E„'5.. 



\{[l + KM] 

ccj a=l 
i=l \i=l / 



n 



,(C7) 



where we have divided by the common factor 
Y\l=i (p exp(^j7iQ)^ on both sides. The normal- 



ization N is defined as 
p 



N 



/ Ildy.d^a.y,] W{y)Y 



(C8) 



This identity holds if and only if the coefficients of lin- 
early independent products of factors of Ki{^a) are iden- 
tical on both sides. Since is normalized, the iden- 
tity holds for Ki{fia) = 0, i.e., /^a = 0. We now 
consider products over factors Ki with i ranging over 
1 < i < p — 1, since Kp — — Y^^iZi ^i- Consider a prod- 
uct with ki factors Ki{^a) (with all different). The 
coefficient on the left-hand side is obtained from combina- 
toric considerations: A factor of Ki either comes directly 
from a term (1 4- Ki) in (C7|, or it results from a term 
[1 + Kp), upon replacing Kp — — Y^^iZi Ki- There are 
k \ 

* j — ki\/ri\{ki — ri)\ different ways to have factors 
ri J 

of the latter origin (each contributing a factor (—1) to 
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the coefficient) and ki — Vi of the former. Then, ki of the indices i to unconstrained rephca a 
(nai) /^-indices with ia = i are already assigned, while 
the remaining na^ — ki indices i need still to be assigned 

to a subset of the n — J^^Z^ ki replica with yet unfixed ia ■ YVi=i {nai ) ! 

The number of possibilities to make disjoint assignments 
for all indices i — 1, ■■■,p is 



(CIO) 



Putting all elements together, the sought coefficient fol- 
lows as 



{nap - ^f^i ny. riLi (""i -h + ny. 



fcl 

(C9) y: 



ri=0 Tp 



^p(nap - X;Li riLi ('^"^ -h + r,)! 



This is normalized by the number of assignments of nai 



On the other hand, the coefficient on the right-hand side is given by 



.Y{U{narr~' 
n\ 



n(-ir 



1=1 



(Cll) 



•' 1=1 \i=l ' —1 



A,;(y)- Ap(y) 



W{v)\ 



ki (S^P „yi\"~J2Ul ki 



nr=i dy. ^(sr=i o^^y^) nr=i (e^- - e^-^)^-- (eLi 

rc«nr=id2/.<5(ELia.2/.)(ELie^-)^ 



(C12) 



Our task is to show that C{fe.} = ^{fe }■ note that a priori Cj^j is only defined for integer and positive nai, while 
^{k } only defined for n < 0, but not necessarily integer. We will show that C'^f.,^ has an analytic continuation 
to positive n and nai which indeed coincides with C^ki} where the latter is defined. Thus we interpret C'j^ j as the 
analytical continuation of the replica expression, which can then be continued to n f 0. 

Let us proceed by computing the numerator in Eq. (C12| (recalling that everywhere we assume X^iLi '^i ~ 1) 



p , p-i 



/oo y y y ^ / y 

\{dy,5[Y,a,yA\{{ey^-ey^f^[Y,. 
-°Oi=l i=l i=l \j=l 

= / n dy.' dypSiY. ^^y^ + yp) e"'" 11 (^'' - 1)'' U + E 

J-OO j^;^ ^ i^^ \ 

poo p-1 p-1 

= / n^y.' n e-""'^^(e^'-l)* 



i=l 

^ p-1 \ 



i=l i=l 
1 



dA 



'i=i 
p-1 



Now we change variables to Oj — and expand the powers 

1 



n dyl e-Hi+^?- ^"0 'Jl [e--.^; (e^^ - 1)^ 

'i=l i=l 



(013) 



{k^} 



ri=0 



-ri— nai — Ig — Aai 



1 j"^ dXe ^ ji f \ i~lV' r(fc, -ri- nai) 

ECi ki) Jo Ai+«-5:&' l\ V * / A'^-'^-"". 



r(-n 



ri=0 rp_"llo r(-n + ELl ^i) i=l 

Finally, we use the relation T{x) — sin(jra;)r(i-z) rewrite this (using that ki and are integers) as 

^p-i 7, ^ P-1 



(014) 



{ki} T-rp 



(-l)P-isin(n7r)/7r 



- V — 

nLi sin(na,^)/^ ^^^^^^^ r(l + na^ - Ef=i n) r(l + na, - A:, + r,) V 



m+n-EUh) 



n (-1)'- • (CIS) 
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The ratio of F functions in Eq. (C15) can be continued 



to positive n. When n and all nai become integers, the 
latter can be written as 



(^-ESfcQ! 

{nap - X^fJi nfr^ {nai - h + r^)l 



(C16) 



Dividing by the normalization factor yields indeed Cjfc.j., 
which completes the proof. 



Note that the normalization factor in the denominator 



of Eq. (C12 1, for n — >■ is given by 



N{n^0)^^^f:^p^[l + O{n)], (C17) 
lli=i sm{naiTT)/TT 

which tends to ^ l/[{-n)P^^ Di when n f 0, as 
calculated previously in Eq. ( 26 ) . 
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